Ab initio study of magnetic structure and chemical reactivity of Cr 2 03 

and its (0001) surface 
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We present the first ab initio density functional theory study of the oxygen-terminated Cr203 (0001) surface 
within the local spin-density approximation (LSD A). We find that spin plays a critical role for even the most 
basic properties of Cr2 03 such as the structure and mechanical response of the bulk material. The surface 
exhibits strong relaxations and changes in electronic and magnetic structure with important implications for 
the chemical reactivity and unusual spin-dependent catalytic activity of the surface. Unlike the bulk, the 
outermost chromium bilayer is ferromagnetically ordered, and the surface oxygen layer exhibits appreciable 
net spin polarization in the opposite sense. Surprisingly, despite this ferrimagnetic order, the chemically 
important states near the Fermi level exhibit ferromagnetic order and thus favor electronic spin alignment of 
species interacting with the surface. Finally, we also find a high density of unoccupied electronic surface states 
available to participate in the chemical reactivity of the surface. 
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I. INTRODUCTION 

Chromium oxide is useful as a catalyst in many ap- 
plications including internal combustion engine emission 
technology JjJ and the dehydrogenation of alkanes [^],[| . 
It is also the primary constituent of passive films pro- 
tecting stainless steels and other high performance in- 
dustrial alloys ||-||. Moreover, the surface is used in 
industry to catalyze the conversion of para- to ortho- 
hydrogen 0| . Intriguingly, this conversion has been ob- 
served to decrease by nearly an order of magnitude as the 
temperature increases through the Neel magnetic order- 
ing temperature ||. That this conversion is so efficient 
for a non-magnetic material such as &2O3 has remained 
a subject of discussion for several decades jjj. Despite 
these intriguing and important phenomena, the surface 
electronic structure of Cr 2 03 is the least studied among 
the transition metal oxides || . 

Studies of corundum-structure metal oxides to date 
have neglected the electronic and magnetic structure of 
the surfaces of Cr203. Manassidis et al. |l(J were one 
of the first to use ab initio density functional theory 
(DFT) within the local density approximation (LDA) to 
study the structure and energetics of basal plane sur- 
faces of a corundum crystal (AI2O3). They found very 
large surface relaxations which resulted in a reduction 
of the surface energy by a factor of two, but, as AI2O3 
is non-magnetic, whether transition metal systems ex- 
hibit this behavior remained unexplored. Veliah et al. 
p[ employed DFT with both the local spin-density ap- 
proximation (LSDA) and non-local LSDA (NLSDA) to 
study Cr m O„ clusters, but only up to m = 2,n = 3, 
too small to provide insight into the behavior of surfaces. 
Classical molecular dynamics calculations p2[ as well as 



ab initio periodic unrestricted Hartree-Fock (UHF) cal- 
culations of bulk Cr203 jl3| and its surface |L4j have been 
carried out, but these studies have neglected the oxygen- 
terminated surface, which is important in oxygen-rich en- 
vironments. Such surfaces have been studied, but in the 
related transition-metal material cv-Fe203. Using spin- 
DFT within the generalized gradient approximation and 
a full potential linearized augmented plane wave (FP- 
LAPW) basis set, Wang et al. jllj studied the a-Fe 2 3 
(0001) surface as a function of increasing oxygen pressure 
and demonstrated that, indeed, the oxygen-terminated 
surface eventually becomes the most stable. 

To better understand the unique properties of Ci^Os, 
we explore the oxygen-terminated (0001) surface of 
Cr203 using the ab initio plane- wave pseudopotential 
approach within LSDA. Section briefly reviews our 



methodology, and Section [III A| gauges the reliability of 
this approach through detailed comparisons with exper- 
imental information on the lattice constant, bulk modu- 
lus, and atomic arrangement and ma gnetic structure of 
the bulk material. Finally, in Section III B , we contrast 
the atomic, electronic, and magnetic structures of the 
surface with those of the bulk. 



II. METHODOLOGY 

Our calculations are carried out within the ab ini- 
tio pseudopotential density-functional formalism which 
has been applied successfully in the past to a wide va- 
riety of solid state and surface systems. (See Q for 
a review.) The present calculations involve three pri- 
mary approximations: the local density approximation to 
density functional theory, the pseudopotential ("frozen- 



1 



core") approximation and the supercell approximation. 

To evaluate the efficacy of local density functional 
in treating the present material, we have calculated the 
properties of bulk 0^03 within both LDA and LSDA 
as parameterized in |n| and jl8| , respectively. These re- 
sults (Section [II A ) show that LSDA gives a description 
far superior to the LDA and also gives results better than 
those reported in the literature for unrestricted Hartree- 
Fock (UHF) calculations |]l3). Accordingly, we employ 
LSDA exclusively for the surface calculations below. 

To represent the chromium and oxygen ionic cores in 
this study, we employ pseudopotentials of the non-local 
Kleinmann-Bylander separable form as generated us- 
ing the optimized pseudopotential procedure of Rappe 
et al. J2(]] . Preliminary calculations carried out on small 
atomic clusters indicate that the states of the argon shell 
are relatively close in energy to the 2s states of oxygen. 
Out of concern for the possible mixing of the argon-shell 
electrons of chromium with the oxygen valence shell, we 
include the argon-shell electrons explicitly in our calcula- 
tion as valence electrons and pseudize only the neon core 
of chromium in our pseudopotential. Convergence tests 
show that both this pseudopotential and that for oxygen 
are well- converged at a plane wave cutoff of 80 Rydbergs, 
which we use for all results reported below. We find that, 
relative to the argon-core pseudopotential, the neon-core 
potential reduces errors in the band structure relative to 
the results of LAPW calculations by a factor of two. 

For the bulk calculations described below, we employ 
the ten-atom corundum primitive cell of Cr 2 03. For this 
cell, we sample the Brillouin zone with an eight k-point 
Monkhorst-Pack |2^| mesh, which folds to four points 
under time-reversal symmetry. This corresponds to a 
reciprocal-space sampling density of 0.37 bohr -1 , which 
is commonly used to treat silicon and which we expect 
to be more than sufficient for this even more ionically 
bonded system. 

For the surface calculations, we employ the supercell 
approach with periodic boundary conditions in all three 
dimensions, thereby allowing us to study isolated slabs 
of Cr203 material with two (0001) surfaces separated by 
7 A of vacuum. This supercell contains five bilayers of 
chromium and six layers of oxygen. We arrange these 
to form a slab with two identical oxygen-terminated sur- 
faces in a unit cell with the in-plane lattice constant held 
at the ab initio value. Such slabs were found sufficient 
for studies of the corresponding surface of Fe2 03 [EBJ and 
avoid the formation of a net electric dipole in the super- 
cell. For the surface cells, we use four k-points (which 
fold to two) in the plane perpendicular to the c-axis of 
the cell. This yields a reciprocal space sampling compa- 
rable to that of our bulk calculations. 

Within the above approximations, we determine the 
quantum state of the system by minimizing the total en- 
ergy over all possible sets of orthonormal electronic wave- 
functions using the analytically-continued functional ap- 
proach 16 2^), which has been recently shown to be 
nearly optimal in a rigorous mathematical sense fl24] . The 



calculations are performed within the DFT++ formal- 
ism |25| , |2"6|l which makes different physical descriptions 
(LDA, LSDA, SIC, Hartree-Fock, etc.) relatively sim- 
ple to explore and provides both high portability and 
good computational performance. Finally, when deter- 
mining structures, we relax the atomic coordinates until 
the Hellmann-Feynman forces on the atomic cores are 
less than 0.03 eV/A in each coordinate direction. 



III. RESULTS AND DISCUSSION 



A. Bulk Cr 2 3 




FIG. 1. Schematic projection of the metallic atoms in the 
corundum structure, viewed atop the (0001) surface, showing 
chromium atoms only. The chromium atoms in the structure 
are colored according depth and are arranged in a sequence 
of bilayers with honeycomb structure (three are shown) . The 
upper and lower chromium atoms of each bilayer appear as 
large and small circles, respectively. 



The crystal structure of Cr2 03 is important for un- 
derstanding the relaxation of the internal coordinates of 
the unit cell, the magnetic order of the bulk, and the 
structure of the surface. Cr203 assumes a corundum- 
type structure (space group R3 c), which has a ten-atom 
primitive cell that is equivalent to a thirty-atom cell on 
a hexagonal lattice. The structure consists of alternating 
oxygen layers and chromium bilayers stacked along the o 
axis of the hexagonal lattice. Figure [j] shows the arrange- 
ment of the chromium bilayers. Each bilayer consists 
of two perfectly planar triangular lattices arranged so 
that their combined projection forms an ideal honeycomb 
structure whose two-element basis consists of one atom 
from each sublayer. The difference between the crystal- 
lographic parameter "z(Cr)" j27j and one-third measures 
the separation of the planes which make up the bilayer. 
In the structure, the chromium bilayers align in an ABC 
stacking sequence along the c-axis such that atoms from 
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the furthest sublayers of each bilayer pair align directly 
along the c-axis. All oxygen atoms in a layer lie in ex- 
actly the same plane and form an approximate triangu- 
lar lattice. The difference between the crystallographic 
parameter "x(O)" p7| an d one-third measures a slight 
contraction of those triangles in this lattice which are 
centered on c-axis aligned pairs of chromium atoms. 

It is known experimentally that, unlike Fe203 and 
V2O3, both AI2O3 and Cr 2 03 contract nearly isotropi- 
cally (within tenths of a percent) under hydrostatic pres- 
sure p3j. Accordingly, for our LDA and LSDA calcula- 
tions, we determined bulk moduli and unit cell volume 
by holding the ratio c/a between the hexagonal crystal 
axes constrained to the experimental value of 2.7407 |||, 
p. 469]. At each volume, we optimized the internal co- 
ordinates of the unit cell, and applied a standard Pulay 
stress correction f|c],|3l| for changes in finite-basis size to 
the resulting total energy. The minimum energy, location 
of the minimum and the second derivative of the energy 
as a function of lattice constant give the cohesive energy, 
unit-cell volume and bulk modulus, respectively. 




(a) (b) 

FIG. 2. LSDA predictions for a thirty atom cell of bulk 
Cr2C>3. Left panel (a) shows a single isosurface of the pre- 
dicted total electron density (n^ + n\), with the surface col- 
ored according to proximity to chromium atoms (white) and 
oxygen atoms (charcoal). Right panel (b) shows the spin po- 
larization, n-f — nj_, for bulk Cr2C>3 within LSDA (thirty atom 
calculation): contour of net up-spin (white), contour of net 
down-spin (charcoal). For both panels some atoms are cut-off 
at the cell boundary, revealing darker shaded interiors. 

Table I compares our results for bulk modulus, unit cell 
volume, cohesive energy (atomization energy to neutral 
atoms), and internal cell coordinates with experiment 
and unrestricted Hartree-Fock results ]l3| . Whereas the 
LDA calculations are in error by about -10% and 44% 
for unit cell volume and bulk modulus, respectively, the 
corresponding errors for UHF are +5% and +13% and 
for LSDA are -1.7% and within experimental error (5%), 
respectively. Moreover, in going from LDA to LSDA, the 
error of the cohesive energy, determined from Born-Haber 
cycle data for Cr 2 03 |||] and the bond strength of molec- 



ular oxygen ||, reduces from -12% to -6%. (Ref. @ 
does not provide a value of the cohesive energy.) 

We find that the LDA prediction of the internal coor- 
dinates of the unit cell is in significant error. We note 
that x(O) for the LDA structure is nearly one-third, in- 
dicating a tendency toward a perfect triangular lattice 
in each O layer. Our LSDA results, on the other hand, 
are in much better accord with experiment, showing even 
somewhat better agreement than reported for UHF [jl3| . 
The uncertainties quoted for x(0) and z(Cr) in the ta- 
ble represent small differences in position of equivalent 
atoms due to numerically imperfect relaxation. 

Given the success of the LSDA description and its clear 
superiority over the LDA (and UHF), all results reported 
below are computed within LSDA. 

Next, we turn our attention to the magnetic struc- 
ture of bulk Cr203. Figure |^b presents our prediction 
for the ground-state spin density. Consistent with ex- 
pectations, we find that the spin polarization in the bulk 
material concentrates almost entirely on the chromium 
atoms (large dimpled contours in the figure). The spin 
order is antiferromagnetic among nearest neighbors in 
each chromium bilayer, leading to an intra-layer order 
in which the spin direction is constant in each sublayer. 
Moreover, we find an inter-layer order such that the spin- 
up/spin-down sequence along the c-axis is identical for 
each bilayer. These results are in complete accord with 
the magnetic order inferred from neutron diffraction ex- 
periments |54j. Finally, we find that the oxygen atoms 
have lobes of opposite spin-polarization with negligible 
net spin. 
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FIG. 3. Antiferromagnetic soliton in Cr203 within LSDA 
(same conventions as in Figure ^b). 

To confirm that the state shown in Figure ^> as the 
ground state within LSDA, we relaxed our initial elec- 
tronic state from completely random wavefunctions. This 
resulted in a spin state in which the magnetic order 
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within each bilayer was the same as described above, 
suggesting that the antiferromagnetic intra-layer order 
is quite stable. However, in this calculation, the spin- 
up/spin-down sequence of the bilayers reversed from one 
ground state to the spin-flip related phase. (See Fig- 
ure |3[) This corresponds to the presence of two anti- 
ferromagnetic solitons in the supercell. LSDA ascribes 
a positive energy to this excitation, indicating that the 
spin ordering shown in Figure ||d represents the ground 
state. From our results, we extract the first ab initio es- 
timate for the (0001) c-axis antiferromagnetic coupling 
constant, Jab ~ 150 cm . Experiments provide highly 
variable estimates for this coupling constant from about 
Jab « 250 cm- 1 to w 350 cm" 1 M. 
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(b) 

FIG. 4. Local chemical softness map of bulk Cr2 03. Por- 
trayed are the HOMOs (a) and LUMOs (b) as calculated 
within LSDA. The figure shows a single isosurface of local 
softness, which is colored according to proximity to either 
chromium (white) or oxygen (charcoal) atoms. 



Finally, Figure f| shows contour surfaces of the elec- 
tron density of the highest occupied molecular orbitals 
(HOMOs) and the lowest unoccupied molecular orbitals 
(LUMOs), representing the local chemical softness [ |36| . 
Species which tend to accept or donate electrons tend 
toward regions with strong HOMO or LUMO concen- 
tration, respectively. The figure shows that, in bulk, 
the chromium atoms interact with either type of species, 
whereas the oxygen atoms are relatively inert. 



B. Oxygen-terminated (0001) surface 

Comparing the relaxed surface slab to the bulk, we find 
that the primary structural change is a significant motion 
of surface oxygen atoms (layer 0[1] in Figure ||a) toward 
the outermost sublayer of chromium atoms (layer Cr[2a]). 
This motion involves both a vertical component perpen- 
dicular to the surface and a lateral component in the 
surface plane. The vertical component reduces the inter- 
planar spacing of layers 0[1] and Cr[2a] by 33%. The 



lateral component rotates triangles of oxygen atoms cen- 
tered above the chromium atoms of layer Cr[2b] by 6.5°, 
similar to the rotation of 10° reported for Fe203 [ p"5| . 
The second most significant change is an increase of 9%, 
or 0.12 A, in the distance between the chromium layer 
Cr[2b] and the next deeper oxygen layer 0[3]. All shifts 
deeper within the structure are less than 0.06 A. Table || 
summarizes our results. 




0[1] 
Cr [2a] 
Cr [2b] 

0[3] 

Cr [4a] 
Cr [4b] 

OH 
Cr [6a] 
Cr [6b] 

0[7] 
Cr [6a] 
Cr [8b] 

0[9] 
Cr[10a] 
Cr[10b] 
0[11] 




(a) (b) 
FIG. 5. LSDA predictions for a twenty-eight atom surface 
slab of Cr2 03: electron density (a) and spin polarization (b) 
(same conventions as Figure ti.) 



The magnetic structure of the surface is significantly 
different from that of the bulk (Figure |b). The most 
striking feature is that the outermost chromium bilayer 
is now ferromagnetically ordered. In addition, we find a 
noticeable spin moment on the outermost oxygen layer in 
the opposite spin orientation, similar to what has been 
found in Fe203 |ljj. From the figure, it is clear that the 
origin of the net spin of the outer oxygen atoms is the 
de-population of the lobes of opposing spin, which are 
normally filled in the bulk. 

In order to understand the chemical reactivity of the 
Cr 2 03 surface, we compare the HOMO and LUMO states 
of the surface slab in Figure §t and |b. These fi gures 
show a remarkable spatial separation between the HO- 
MOs and LUMOs. The HOMO states, which would in- 
teract most strongly with electron-accepting species such 
as protons, contract into the bulk and leave almost no 
concentration on the outer oxygen atoms. On the other 
hand, the LUMO states, which interact with electron- 
donating species such as chloride or sulfide ions, concen- 
trate strongly on the surface. This suggests that such 
electron-donating species would be more disruptive to 
the surface than would be their electron-accepting coun- 
terparts. 
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FIG. 6. Local chemical softness map for the surface slab of 
Cr 2 3 : HOMOs (a), LUMOs (b). Colors indicate proximity 
to chromium (white) and oxygen (charcoal). 



j* » D * « P I 1 




(a) 



FIG. 7. Spin-dependence of the local chemical softness map 
for the surface slab of Cr 2 3 : HOMOs (a), LUMOs (b). Col- 
ors indicate spin direction. 



hydrogen which has been observed to occur more readily 
in the presence of magnetically ordered surfaces of Cr 2 03 



IV. CONCLUSIONS 

The preceding results for the bulk crystal establish the 
critical role of spin in the physics and chemistry of C^OV 
We show that the plane- wave pseudopotential LSDA ap- 
proach gives a good description of the structure, mechan- 
ical response and magnetic order of the bulk. Through 
this approach, we provide the first ab initio estimate of 
the antiferromagnetic coupling constant. 

The oxygen-terminated (0001) surface of Cr 2 03 ex- 
hibits strong inward relaxation of the outermost oxygen 
layer with an associated rotational reconstruction, as has 
been found in previous theoretical studies of similar ma- 
terials. Understanding surface processes such as catalysis 
and corrosion initiation, however, requires fundamental 
knowledge of the electronic structure and cannot be in- 
ferred from the surface relaxation alone. 

We find that the electronic structure of the surface 
differs radically from that of the bulk. In particular, the 
outermost chromium bilayer does not exhibit the anti- 
ferromagnetic order of the bulk but is ordered ferromag- 
netically. Moreover, the outer oxygen layer, rather than 
exhibiting zero net spin polarization, develops an appre- 
ciable net spin, but in the sense opposite to that of the 
underlying ferromagnetic chromium bilayer. Addition- 
ally, the surface manifests a high density of low-energy 
unoccupied electronic states available to influence chemi- 
cal reactions at both oxygen and chromium surface sites, 
but with negligible penetration into the bulk. There is 
also some density of high-energy occupied states reach- 
ing outward from the bulk into the outermost chromium 
bilayer. Remarkably, although the spins of the outer- 
most oxygen and chromium atoms orient in opposite di- 
rections, the chemically relevant electronic states create a 
blanketing environment of uniform spin which is available 
to encourage spin alignment of adsorbed species. This 
new understanding provides a path for insight into novel 
magneto-chemical effects and their implications for catal- 
ysis. 
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Model B (Mbar) V (A 3 ) E (eV) 1/3 - x(Q) z(Ci) - 1/3 

LDA 3.39(10) 43.34(04) 23.63 -0.004(2) 0.00626(1) 

LSDA 2.35(12) 47.19(14) 25.24 0.03133(8) 0.01647(1) 

UHF 2.66(10) 50.5 — 0.03273 0.01722 

Published Experimental Values 
Ref. M 2.38(12) 47.997(5) = = = 

Ref. [28| 2.31(30) 48.0 0.0282 0.01437 

Ref. [ 32 1 26.87 1 " 

Ref. H 0.02763 0.01417 

' Computed with Born-Haber cycle. (See Text.) 

TABLE I. Comparison of computed bulk moduli B, vol- 
ume per formula unit V, cohesive energy E, and internal cell 
coordinates x and z in ten-atom cell of Cr2 03 with experimen- 
tal values. Volume and cohesive energies (relative to isolated 
atoms) are per Cr2C>3 formula unit. 



Layer pair 


Bulk 


Surface 




spacing [A] 


spacing [A] 


0[1]-Cr[2a] 


0.94 


0.63 


Cr[2a]-Cr[2b] 


0.39 


0.40 


Cr[2b]-0[3] 


0.94 


1.05 


Q[3]-Cr[2b] 


0.94 


0.92 



TABLE II. Ab initio prediction of changes in interlayer 
spacing at the (0001) a-Cr2C>3 surface. 
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